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Abstract 

The  method  of  thermodynamic  simulation  was  used  to  study  chemical  interactions  in  the  cathode  half-cell  of  lithium-ion  battery,  which  contained 
LiCo02  as  cathode  active  material  and  a  LiPF6  solution  as  the  electrolyte.  It  was  shown  that  in  the  temperature  range  298-400  K  in  thermodynamic 
equilibrium  state  a  layer  of  solid  products  of  chemical  reaction  formed  on  the  cathode/electrolyte  interface.  The  layer  predominantly  consisted  of 
LiF  and  LiP03  in  the  molar  ratio  about  2:1  (corresponding  to  the  volume  ratio  2:3).  In  equilibrium  state  also  some  soluble  interaction  products 
formed,  namely:  CoFi  (considerable  quantity).  POF3  and  PF5.  The  concentration  of  the  later  two  substances  was  small,  but  it  increased  with 
increase  of  the  temperature.  Furthermore,  in  equilibrium  state  oxygen  gas  formed  as  a  product  of  chemical  interactions  in  cathode  half-cell. 

Thermodynamic  characteristics  of  C02O3  and  LiCoCT  compounds  were  determined  using  a  set  of  calculation  methods. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  study  of  processes  concerned  with  the  degradation  of 
lithium-ion  batteries  is  one  of  the  main  research  topics  in  the 
field  of  chemical  power  sources.  In  spite  of  the  fact  that  the 
characteristics  of  the  cells  produced  by  different  manufacturers 
vary  within  broad  range,  there  are  several  common  problems  for 
Li-ion  batteries  (hereafter  LIB),  namely:  capacity  fading,  aging, 
acceleration  of  these  processes  with  increase  of  the  temperature, 
and  considerable  gassing  under  heating  [1-5]. 

Various  aspects  of  degradation  phenomena  in  LIB  have  been 
considered  in  a  number  of  papers  and  the  relevant  results  have 
been  summarized  in  reviews  [6-8].  Special  attention  was  paid  to 
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the  analysis  of  electrode  surface  layers.  Most  of  the  studies  deal 
with  the  formation  of  anode  surface  layers  (see,  e.g.  [9-12]), 
while  the  cathode/electrolyte  interface  was  much  less  studied 
[13-15], 

The  present  study  is  dedicated  to  chemical  interactions  in  the 
cathode  half-cell  of  LIB  and  to  the  processes  concerned  with 
these  interactions  and  leading  to  degradation  of  the  cells  in  the 
course  of  long-term  storage.  The  study  was  carried  out  using  the 
cells  produced  by  Samsung  SDI  on  a  mass  scale,  which  contain 
LiCoCE  as  cathode  active  material  and  1 . 1 5  M  LiPFf,  solution 
in  mixture  of  alkyl  carbonates  as  the  electrolyte.  The  Part  I  of 
the  series  of  papers  reports  results  of  a  thermodynamic  simula¬ 
tion  of  the  interactions  on  the  cathode/electrolyte  interface.  The 
Part  II  is  concerned  with  the  experimental  study  of  chemical 
processes  taking  place  in  the  cathode  half-cell.  X-ray  diffrac¬ 
tion  (XRD),  Fourier  transform  infrared  spectroscopy  (FTIR), 
X-ray  photoelectron  spectroscopy  (XPS)  and  electrochemical 
impedance  spectroscopy  (EIS)  were  used. 
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2.  Method  of  thermodynamic  simulation 

Theoretical  principles  of  the  thermodynamic  analysis  of  com¬ 
plex  reagent  systems  were  formulated  in  the  classical  ther¬ 
modynamics  by  Gibbs.  Lewis  and  Randall  [16]  bridged  the 
theory  and  the  practice.  One  of  the  pioneering  algorithms  for 
computer-aided  calculation  of  the  equilibrium  composition  was 
proposed  by  Brinkley  and  Kandiner  [17,18].  Another  algo¬ 
rithm,  which  is  based  on  the  minimization  of  Gibbs  energy, 
was  developed  later  [19].  The  first  software  kit,  which  was 
designed  for  mass  calculations  of  the  equilibrium  composi¬ 
tion  and  had  a  database  of  thermodynamic  properties  of  indi¬ 
vidual  substances,  was  elaborated  by  Zeleznik,  Gordon  and 
McBride  [20,21],  The  methods  of  thermodynamic  simulation 
proved  to  be  incredibly  fruitful  for  the  analysis  of  metallur¬ 
gical  processes.  This  promotes  further  improvement  of  the 
simulation  methods  [22-24],  By  now,  several  hundreds  of 
algorithms  and  software  kits  have  been  developed  to  calcu¬ 
late  the  equilibrium  composition  of  multi-component  systems. 
The  most  popular  algorithms  are  comprehensively  reviewed  in 
[25-28], 

The  method  of  thermodynamic  simulation  (TS)  essentially 
consists  in  a  comprehensive  thermodynamic  analysis  of  the  equi¬ 
librium  state  of  the  system  as  a  whole.  Modern  algorithms  are 
based  on  the  variation  principles  of  the  thermodynamics  devel¬ 
oped  by  Dorn  [29].  It  is  assumed  that  the  thermodynamic  system 
is  closed  and  isolated,  that  the  system  obeys  the  laws  of  statisti¬ 
cal  physics  (it  is  true,  if  the  number  of  atoms  in  the  system  is  over 
105),  that  the  system  is  isotropic  and  not  influenced  by  magnetic, 
electric  and  gravitational  fields.  A  full  thermodynamic  analysis 
of  the  equilibrium  state  of  a  model  system  allows  one  not  to 
take  into  account  individual  chemical  reactions.  The  final  result 
of  the  simulation  is  the  quantitative  determination  of  equilib¬ 
rium  interaction  products  in  a  given  model  system.  Since  the 
aim  of  this  study  is  to  determine  the  degradation  mechanism  of 
LIB  during  a  long-term  storage  one  can  assume  that  the  system 
under  the  study  is  close  to  equilibrium  state.  The  comparison 
between  calculated  and  experimental  data  makes  the  interpreta¬ 
tion  of  experimental  data  much  easier  since  only  the  products 
of  interaction,  which  are  thermodynamically  allowed  and  there¬ 
fore  should  form  according  to  modeling,  could  be  taken  into 
account. 

All  modeling  calculations  were  done  using  ASTRA  4  soft¬ 
ware  kit  with  the  ASTRA.BAS  database  [30].  In  the  frame  of 
the  model  used  the  criterion  of  the  equilibrium  state  was  the 
maximum  of  the  entropy  of  the  system.  The  following  ther¬ 
modynamic  values  of  the  studied  compounds  were  used  for 
the  simulation:  standard  formation  enthalpy  (Af/L?9g),  stan¬ 


dard  entropy  (S^g),  heat  capacity  under  standard  conditions 
(c°o98),  the  incremental  enthalpy  (//|9g  —  Hq),  and  the  coef¬ 
ficients  of  the  temperature  dependence  of  the  heat  capacity 
in  Mayer-Kelli’s  equation  cp  =  a  +  bT—  cT~- .  Thermodynamic 
values  of  some  compounds  could  be  found  in  the  literature 
and  added  to  the  database.  Thermodynamic  values  of  LiPF6 
were  taken  from  [31,32].  Characteristics  of  LiCo02  and  C02O3 
were  calculated  using  empirical  methods  [33]  and  the  method 
of  thermodynamic  simulation  [34]  adopting  the  model  of  ideal 
solutions  of  interaction  products  [35].  Reference  characteris¬ 
tics  necessary  for  calculation  of  thermodynamic  properties  were 
taken  from  [36,37]  and  ASTRA.BAS. 

3.  Results  and  discussion 

3.1.  Calculation  of  thermodynamic  properties  of  the 
components  of  the  system 

3.1.1.  Calculation  of  thermodynamic  properties  of  C02O3 

Thermodynamic  characteristics  of  LiCo02,  which  are  not  yet 

reported  in  the  literature,  could  be  calculated  using  thermody¬ 
namic  properties  of  cobalt  (III)  oxide.  However,  there  exists  no 
reliable  literature  data  concerned  with  C02O3  as  well.  There¬ 
fore,  thermodynamic  values  of  C02O3  were  calculated  by  the 
empirical  method.  The  calculations  were  based  on  the  analogy 
of  thermodynamic  properties  between  cobalt  (CoO,  C02O3  and 
C03O4)  and  iron  (FeO,  Fe203  and  Fe30/t)  oxides.  The  stan¬ 
dard  enthalpy  of  formation  and  standard  entropy  of  Fe304  is 
larger  than  the  corresponding  total  values  of  the  mixture  of 
FeO  +  Fe203  by  22kJmol“l  and  2JK_1mol_1,  respectively 
[34],  Assuming  the  same  relationship  between  these  values  for 
cobalt  oxides,  we  estimated  the  standard  enthalpy  of  formation 
and  standard  entropy  of  C02O3  using  well-known  thermody¬ 
namic  characteristics  of  CoO  and  C03O4. 

It  is  known  that  for  some  transition  metals  (including  iron) 
cp  of  metal  oxide  depends  linearly  versus  index  x  in  the  oxide 
chemical  formula  MeOA  [38].  Assuming  that  this  dependence 
is  correct  for  cobalt  oxides,  we  estimated  the  heat  capacity  of 
C02O3  under  standard  conditions  (dependence  of  the  heat  capac¬ 
ity  versus  index  x  in  the  formula  CoO A  was  linearly  extrapolated 
tojc=  1.5). 

The  enthalpy  increment  and  the  coefficients  in  Mayer-Kelli’s 
equation  (describing  dependence  of  the  heat  capacity  versus 
temperature)  were  calculated  using  empirical  equations  given 
in  [39]. 

Calculated  thermodynamic  values  of  C02O3  are  given  in 
Table  1. 


Table  1 

Calculated  thermodynamic  properties  for  C02O3  and  LiCo02 


Compound 

Af ff2°9g  (kj  mol-1) 

S?98  (JK-1  mor1) 

cp=a  +  bT— 

cT~2  (JK-'mol-1: 

) 

H29S  ~  Ho  (Jmor1) 

a 

bx  103 

c  x  10~5 

C02O3 

-631.6 

62.3 

113.3 

4.5 

46.46 

16000 

LiCo02 

-738 

84.4 

90.64 

0.47 

24.0 

9800 
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3.1.2.  Calculation  of  thermodynamic  properties  ofLiCo02 

The  standard  enthalpy  of  formation  (Af H%9S),  the  standard 
entropy  (S^98)  and  the  heat  capacity  under  standard  condi¬ 
tions  (c°  298 )  of  the  lithium  cobaltite  were  calculated  by  the 
method  proposed  by  Moiseev  and  Vatolin  [34],  The  method  is 
based  on  the  model  of  ideal  solutions  of  interaction  products 
(which  is  a  particular  case  of  the  model  of  ideal  associated  solu¬ 
tions  proposed  by  Prigogine  and  Defay  [40])  and  includes  TS. 
The  method  is  described  in  details  in  [34,41].  This  empirical 
method  for  calculation  of  thermodynamic  characteristics  was 
applied  by  Moiseev  and  Vatolin  to  several  hundreds  of  complex 
oxides  and  was  proved  to  give  a  good  fit  with  experimental  data 
[34], 

To  simulate  the  properties  of  LiCo02  we  set  the  composition 
of  the  initial  system  as  follows:  1  mol  Li  +  1  mol  Co  +  20  mol  CL. 
All  condensed  phases  in  the  modeled  system  were  assumed  to  be 
components  of  an  ideal  solution  of  interaction  products.  In  other 
words,  all  known  oxides  in  the  systems  Co-O  and  Li-0  were 
considered  as  constituents  of  this  solid  solution.  The  ASTRA.4 
software  package  and  the  ASTRA.BAS  database  were  used  for 
the  thermodynamic  simulation  of  the  equilibrium  composition 
of  the  system  in  the  temperature  range  from  273  to  1500  K  (with 
a  step  equal  to  10  K)  under  total  pressure  105  Pa.  For  each  given 
temperature  T  numbers  of  moles  m;  of  simple  oxides,  which 
form  solution  of  interaction  products,  were  calculated.  The  val¬ 
ues  obtained  were  used  to  calculate  the  oxygen  index  y  of  the 
solution: 

y  =  y^mjXj,  (1) 

i 

where  Xj  is  the  oxygen  index  of  the  /th  simple  oxide.  The  tem¬ 
perature  dependence  of  the  total  oxygen  index  y  for  the  system 
“1  mol  Li  +  1  mol  Co  +  20  mol  O2”  is  shown  in  Fig.  1 .  The  vari¬ 
ation  of  the  total  oxygen  index  y  takes  place  thanks  to  change 
of  the  ratio  between  different  lithium  and  cobalt  oxides  forming 
the  solution  (LLO,  LLCL,  CoO,  C02O3  and  C0O2).  As  could 
be  seen  in  Fig.  1 ,  total  oxygen  index  y  of  the  solution  of  inter¬ 
action  products  coincides  with  oxygen  index  of  LiCoCL  ( y  =  2) 


Fig.  1 .  Temperature  dependence  of  total  oxygen  index  y  of  an  ideal  solution  of 
interaction  products  for  “1  mol  Li  +  1  mol  Co  +  20  mol  O2”  model  system. 


at  473  K.  The  numbers  of  moles  /«,  of  the  simple  oxides  form¬ 
ing  the  solution  of  interaction  products  at  473  K  were  used  for 
the  calculation  of  thermodynamic  characteristics  of  LiCo02  as 
shown  below. 

The  standard  entropy  of  LiCo02  was  calculated  using  the 
equation  given  in  [34]: 

52°98(LiCo02)  (IK”1  mor1)  =  ]Tm,S2°98(0,  (2) 

i 

where  m,  is  the  number  of  moles  of  /th  simple  oxide  forming 
the  solution  of  interaction  products  and  S%9s(i)  is  the  standard 
entropy  of  the  /th  simple  oxide. 

The  heat  capacity  under  standard  conditions  was  calculated 
using  the  equation  given  in  [34] : 

cp298(LiCo02)  (J  K_1  mol”1)  =  ^m{c°298(0,  (3) 

i 

where  c°pl98(i)  is  the  heat  capacity  under  standard  conditions  of 
the  /th  simple  oxide. 

The  standard  enthalpy  of  formation  of  LiCo02  was  calculated 
using  the  equation  given  in  [34]: 

Af//2°98(LiCo02)(Jmor1) 

=  ^m;Afff298(0  +  ^m,-(min)Afff298(0,  (4) 

i  i 

where  w,- (min)  is  the  number  of  moles  of  /th  simple  oxide,  except 
Li20  and  Co2C>3  (which  are  basic  for  LiCo02)  and  Af//|98(z) 
is  the  standard  enthalpy  of  formation  of  the  /th  simple  oxide. 

The  incremental  enthalpy  (/7298  —  //[,)  and  the  coefficients 
of  the  temperature  dependence  of  heat  capacity  in  the  equation 
cp=a  +  bT  —  cT~2  were  calculated  using  empirical  equations 
given  in  [39]. 

The  calculated  thermodynamic  values  of  LiCoCL  are  given 
in  Table  1 . 

3.2.  The  reaction  zone  model 

The  interaction  between  the  cathode  active  material  (LiCoCb) 
and  the  electrolyte,  which  consists  of  the  mixture  of  alkyl  car¬ 
bonates  and  LiPFg  lithium  salt,  was  simulated  using  the  fol¬ 
lowing  model  of  the  reaction  zone.  It  was  assumed  that  alkyl 
carbonates  were  chemically  inert  to  the  lithium  cobaltite.  It  is 
known  that  all  these  compounds  have  a  wide  electrochemical 
window  (>4.5  V)  and  they  are  stable  with  respect  to  cathode 
materials  [13].  It  was  assumed  therefore  that  all  chemical  pro¬ 
cesses  in  the  cathode  half-cell  are  concerned  with  interaction 
between  LiCoC>2  and  LiPFf,  only.  The  reaction  zone  of  the 
model  system  is  drawn  schematically  in  Fig.  2.  Also  it  was 
assumed  that  in  the  beginning  of  the  chemical  reaction  there 
was  a  large  excess  of  lithium  cobaltite  in  the  system  and  then 
the  ratio  between  the  components  changed  smoothly  and  Li  Ply, 
part  increased  until  the  entire  LiCoC^/electrolyte  interface  got 
covered  with  a  monomolecular  layer  of  the  lithium  salt.  To  cal¬ 
culate  the  composition  of  the  model  system  quantitatively  the 
LiCoC^LiPfy,  ratio  was  changed  by  2-8  mol%  steps.  Thus,  the 
first  step  of  the  simulation  corresponded  to  the  individual  lithium 
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Fig.  2.  Scheme  of  the  reaction  zone  on  the  “LiCo02/electrolyte  solution”  interface. 


cobaltite,  while  the  last  step  corresponded  to  the  interaction 
between  a  monomolecular  Li  Ply,  layer  with  a  monomolecu- 
lar  LiCoCF  layer  counting  per  surface  unit  of  the  interface 
(Fig.  2). 

Thickness  of  the  monomolecular  layer  of  the  ;th  reagent,  <5(i), 
was  calculated  using  the  equitation  given  in  [42,43]: 


<5(i)  (cm)  =  2  R(i)T 


VimjT) 
2: tA 


1  1/3 


(5) 


where  R(i)r  is  the  effective  radius  of  the  interaction  (the  Onsager 
radius)  between  molecules  of  the  zth  reagent  at  a  temperature  T\ 
Vjm(T)  the  molar  volume  of  the  zth  reagent  at  a  temperature  T\ 
and  A  is  the  Avogadro  constant. 

Hence,  the  mass  of  the  initial  reagents  LiCo02  and  LiPFf,  per 
unit  area  of  the  reaction  zone  was  determined  using  the  following 
equation: 


/«(0(gcm  2)  =  2d(i) 


VimjT)' 

2nA 


1/3 


(6) 


where  d(i)  is  the  density  of  the  zth  reagent.  The  calculated  molar 
fraction  of  LiCoCF  at  the  last  step  of  the  simulation  amounted 


0.56. 


3.3.  Thermodynamic  simulation  of  chemical  interactions  in 
the  cathode  half-cell 

To  perform  thermodynamic  simulation,  we  set  the  compo¬ 
sition  of  the  initial  system  as  follows:  LiCoCF  +  LiPFg  +  Ar. 
LiCo02:LiPF6  ratio  was  varied  in  the  range  0-0.56  mentioned 
above.  The  argon  concentration  in  the  model  system  was  fixed  to 
1  mass%.  The  simulation  was  done  for  total  pressure  of  105  Pa 
and  temperatures  range  298-400  K  (with  10  K  steps).  Based  on 
the  preliminary  estimations  we  assumed  that  traces  of  water  in 
the  electrolyte  (50  ppm)  made  negligible  influence  on  the  equi¬ 
librium  state  of  the  system  and  therefore  they  were  not  taken 
into  account. 

The  results  of  the  simulation  for  the  temperature  298  K  are 
shown  in  Fig.  3.  The  increase  of  the  temperature  slightly  effected 
equilibrium  concentrations.  The  main  effect  was  exclusively  due 
to  the  decomposition  of  LiPFg  at  348  K.  In  the  temperature  range 
273-348  K  the  concentrations  of  PF5  and  POF3  increased  from 
10-7  to  10-5  and  from  10-5  to  10-3,  respectively  (mole  fraction 
units). 

The  results  obtained  let  us  conclude  that  solid  LiF  and  C02O3 
as  well  as  small  quantities  of  gaseous  POF3  and  PF5  products 
formed  in  equilibrium  state  on  initial  stage  of  the  interaction 
(which  corresponds  to  the  excess  of  lithium  cobaltite).  The 
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Fig.  3.  The  dependence  of  equilibrium  concentrations  of  interaction  products 
in  “LiPF6  +  LiCo02  +  Ar”  model  system  at  298  K  vs.  the  ratio  of  the  reagents  in 
the  reaction  zone  (in  mole  fraction  units). 

increase  of  LiPF6  concentration  led  to  the  formation  of  C0F2, 
LiP03  and  oxygen  gas  instead  of  the  cobalt  (III)  oxide. 

Let  us  analyze  the  simulation  results.  In  thermodynamic  equi¬ 
librium  state  the  layer  of  solid  interaction  products  formed  on 
the  cathode/electrolyte  interface.  The  layer  consisted  predomi¬ 
nantly  of  LiF  and  LiPO,  in  molar  ratio  about  2: 1  (it  corresponds 
to  volume  ratio  2:3).  One  can  expect,  that  this  layer  can  separate 
the  reagents  and  inhibit  further  chemical  interaction,  but  does 
not  prevent  the  transport  of  lithium  cations  (LiP03  is  a  known 
lithium-conducting  solid  electrolyte  [44]).  However,  along  with 
passivating  film  on  the  surface  of  lithium  cobaltite  particles  a  sol¬ 
uble  compound,  namely  C0F2,  forms.  C0F2  probably  dissolves 
in  electrolyte  solution.  Moreover,  Lewis  acids  POF3  and  PF5 
form.  The  presence  of  these  compounds  in  electrolyte  solution 
even  in  small  quantities  can  cause  the  decomposition  of  organic 
solvents  and  the  formation  of  several  liquid  and  gaseous  prod¬ 
ucts  [45].  Finally  in  cathode  half-cell  gaseous  oxygen  forms, 
which  can  participate  in  the  chemical  interactions  within  the 
anode  half-cell. 

3.4.  Comparison  of  calculation  and  experimental  results 

The  results  of  the  experimental  study  of  the  chemical  pro¬ 
cesses  taking  place  in  the  cathode  half-cell  of  LIB  containing 
LiCoCL  as  cathode  active  material  and  LiPFg  solution  as  elec¬ 
trolyte  are  describe  in  details  in  Part  II  of  this  series  of  the 
papers.  To  compare  TS  and  experimental  results  briefly,  sum¬ 
mary  of  XPS  data  is  given  below.  As  it  was  shown  by  the  study 
of  LiCoCL  surface  after  contact  with  Li  PFf,  solution  the  layer  of 
the  interaction  products  contained  LiF,  LiPCL  and  C0F2.  So, 
experimental  data  are  in  good  agreement  with  the  results  of 
thermodynamic  simulation  of  chemical  interactions  in  cathode 
half-cell  of  lithium-ion  battery. 

4.  Conclusion 

The  study  of  chemical  interactions  in  the  cathode  half-cell  of 
LIB  by  the  method  of  thermodynamic  simulation  demonstrated 


that  LiCoC>2  and  LiPFf,  are  not  chemically  inert  to  each  other. 
The  data  obtained  lets  one  conclude  that  their  interaction  leads  to 
the  formation  of  the  layer  of  interaction  products  on  the  surface 
of  LiCoCL  particles.  The  layer  predominantly  consists  of  LiF 
and  LiPCL.  Considerable  quantity  of  soluble  C0F2  and  small 
amounts  of  Lewis  acids  POF3  and  PF5  also  form  and  dissolve  in 
electrolyte  solution.  The  gaseous  oxygen  can  form  in  the  cell  as 
well.  The  increase  of  the  temperature  slightly  effects  the  equi¬ 
librium  composition  of  the  interaction  products,  except  POF3 
and  PF5,  whose  concentration  rises  considerably  with  increase 
of  the  temperature,  but  still  remains  not  significant. 
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